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ABSTRACT _ _ 

The reconstruction algorithm introduced by |Eisenstein et aL ( 2007| ), which is widely used 
in clustering analysis, is based on the inference of the first order Lagrangian displacement 
field from the Gaussian smoothed galaxy density field in redshift space. The 2smoothing 
scale applied to the density field affects the inferred displacement field that is used to move 
the galaxies, and partially 2erases the nonlinear evolution of the density field. In this arti¬ 
cle, we explore this crucial step 2in the reconstruction algorithm. We study the performance 
of the reconstruction technique using two metrics: first, we study the performance using the 
anisotropic clustering, extending previous studies focused on isotropic clustering; second, we 
study its effect on the displacement field. We find that smoothing has a strong effect in the 
quadrupole of the correlation function and affects the accuracy and precision 2with which we 
can measure Da{z) and H{z). We find that the optimal smoothing scale to use in the re¬ 
construction algorithm applied to BOSS-CMASS is between 5-10 h~^Mpc. Varying from the 
“usual” 15/i“^Mpc to bh~^Mpc 2shows ^ 0.3% variations in Da{z) and ^ 0.4% H{z) and 
uncertainties are also reduced by 40% and 30% respectively. We also find that the accuracy 
of velocity field reconstruction depends strongly on the smoothing scale used for the density 
field. We measure the bias and uncertainties associated with different choices of smoothing 
length. 


1 INTRODUCTION 

The Baryon Acoustic Oscillations (BAO) is without any doubt a 
robust and promising probe for decrypting dark energy (DE). Fur¬ 
thermore, BAO plays a central role in current and future DE exper¬ 
iments devoted to understand cosmic expansion ( |Levi et al.|2013[ 
[Laureijs et al.||20lT| |Spergel et al.||2013^ . The BAO signature 
corresponds to the imprint left by sound waves generated during 
the early Universe in the baryon-photon fluid and propagated until 
the time of decoupling, when the photon decouples from baryons 
and the sound waves get frozen. While the linear physics of these 
baryon acoustic oscillations is well understood, the nonlinear evo¬ 
lution of the matter density field leads to a coupling of the Fourier- 
modes, which generates a damping of the Baryonic Acoustic Oscil¬ 
lations (BAO) in the power spectrum of galaxies as well as a shift in 
the BAO peak position. The damping in the power spectrum trans¬ 
lates to a blurring of the baryonic acoustic peak in the correlation 
function of galaxies in configuration space. This reduction of the 
contrast in the baryonic acoustic feature increases the uncertainty 
in the BAO distance derived from measurements. 

The major contribution to this damping comes from the bulk 
flows that shift the positions of the galaxies from the linear the¬ 
ory prediction ( [Eisenstein et al.|2007| l. Early work on reconstruc¬ 


tion ( [Eisenstein, Seo & White|200^|Padmanabhan, White & Cohn| 

|2009| l suggested that this effect could be partially reversed using 
the density held to infer the gravitational potential that sources the 
movement of the galaxies. Running backwards the gravitational 
flow restores the BAO feature and reduce the errors in the distance 
measurements. 

Padmanabhan et ak] \2012) extended the methodology from 


[Eisenstein et al. { |2007| l by applying reconstruction to galaxy cata¬ 
logues. Since then, reconstruction has been successfully applied to 
galaxy clustering analysis in different galaxy catalogues and sur¬ 
veys. In most of the caseQ reconstruction is shown to increase 
the precision of the measurements in the samples analyse d (jPad- 


manabhan et al.|20T^ [Anderson et al.|2012||20T3]|2014[|Tojeiro 


et al.|2014[ Kazin et al.|2014[[Ross et al.|2014[ l with up to 45 per¬ 

cent improvement in the error bars in certain samples. These results 
are responsible for reconstruction becoming an essential part of the 
clustering analysis. 


^ Few samples in SDSS galaxy catalogues, such as DR9 CMASS sample 
jAnderson et al.|20f2) and the DR 10 LOWZ sample (Tojeiro et al.|2014h 
reported almost no improvement, but such results are consistent with the ex¬ 
pected results using mock galaxy catalogues given the already small initial 
error compared to the average of the pre-reconstruction catalogues. 
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While reconstruction has shown to be a successful mechanism 
to obtain more precise BAO measurements, the performance of re¬ 
construction algorithm is still under intense investigation. Given 
the current precision in BAO distance measurements, the study 
of the associated uncertainties becomes essential for current and 
future surveys. Some efforts to understand reconstruction analyti¬ 
cally have appeared in the literature: [Padmanabhan, White & Cohn| 
( |2009| ) provided an analytic formalism within the context of La- 
grangian perturbation theory; |Noh, White & Padm anabhan| |2009| ) 
extended this formalism for biased tracers. |White( ( [2015| ) developed 
the formalism to describe post-reconstruction correlation functions 
within the Zeldovich Approximation. |Xu et al. \20\2\ stud ied the 


isotropic effects of reconstruction while [Anderson et al.| ( |2012| ) 
studied the anisotropic effects. Reconstruction performance also 
has been studied using N-body simulations ( |Seo et al.|2006|[20T0] ). 
[Padmanabhan et ah] ( |2012| ) provides a first exploration of the ro¬ 
bustness of reconstruction methodology against different imple¬ 
mentation choices in the isotropic BAO analysis. [Burden et al.j 
( [2014[ ) and [Burden et al[ ( [2015[ ) provided a more detailed empiri¬ 
cal study of the systematics (density dependence, geometry effects, 
RSD corrections) in the reconstruction algorithm. Both of these 
studies concentrated on the reconstruction effects on isotropic clus¬ 
tering. [Vargas-Magana et ^ ( [2014[ ) also studied the results post¬ 
reconstruction but focused on anisotropic fitting systematics. 

The reconstruction algorithm requires an estimate of the den¬ 
sity field in order to estimate the displacement field using the Zel- 


dovich approximation ([Zeldovich|1970 

1 . The density field is usu- 

ally smoothed with a Gaussian kernel ( 

Padmanabhan et al.|2012| 


will source the displacement field. A large smoothing scale could 
erase cosmological information, reducing the effect of reconstruc¬ 
tion. Furthermore, information in the regions sourcing the non¬ 
linear growth will be suppressed. On the other hand, a too small 
smoothing scale increases the noise in the linear density field. As 
it is implemented now, the smoothing scale is a free parameter. In 
Table we summarise different studies related to the smoothing 
scale; we show the references, scales explored, the method of anal¬ 
ysis (i.e configuration of Fourier space), and the kind of mocks used 
for the study. [ White| ( [20 1 0| studied the shot noise effect on recon¬ 
struction and found that for low density tracers, a large smoothing 
scale performs better in terms of isotropic clustering as it gener¬ 
ates a smaller shift in the BAO measurement. Most BAO clustering 
analyses ( Eisenstein et al. 12007) [Padmanabhan et al.|2012[ [Ander-[ 


[son et al. 2012 [2013[[2014[ l use a Gaussian smoothing kernel of 

R = 10-20 h~^ Mpc. Some small deviations from this smoothing 
length have been shown not to alter the results (see Appendix B of 
Anderson et al. 2012 and Padmanabhan et al. 2012a). [Burden et aT] 
( [2014[ ) studied the impact of the smoothing length on the isotropic 
BAO analysis. They found that the bias in the measurement of the 
isotropic dilation parameter a is reduced using a smoothing scale of 
Mpc for CM ASS and a smoothing scale of 10/i“^ Mpc for 
LOWZ. In this paper, we extend those analyses to present a study 
of the effect of smoothing scale on reconstruction performance in 
the anisotropic BAO analysis. 

The motivation for this study is illustrated in Figure We 
show the mean and root mean square (RMS) of the correlation 
function from 100 mock catalogues before and after reconstruc¬ 
tion (hereafter, pre- and post-reconstruction). The different colours 
are different kinds of mock catalogues. Darker shades are pre¬ 
reconstruction, and lighter shades are post-reconstruction cata¬ 
logues. In blue are the PTHALOS mock catalogues used for the 
analysis of the Baryonic Oscillations Spectroscopic Survey (BOSS) 


(Dawson et al|2013> o f Sloan Digital Sky Survey III (SDSS-III) 
( Eisenstein et al||2011^ galaxy samples Data Releases 9, 10 and 
11 (DR9, DRIO, DRll). PTHALOS mocks are based on second 
order Lagrangian Perturbation Theory ( [Manera et al.[[2013] l. In 
green are the Quick Particle Mesh mock catalogues used to anal¬ 
yse BOSS DR12 based on low resolution N-body simulations com¬ 
bined with HOD models for populating halos with galaxies. Fi¬ 
nally, in red are the MD-PATCHY (Kitaura et al, 2015; companion 
paper) mocks based on Augmented Lagrangian Perturbation The¬ 
ory and stochastic bias prescriptions generated for BOSS DR 12 
analysis. In all types of simulations, reconstruction succeeded in 
1) sharpening the BAO feature in the monopole, and 2) reducing 
the quadrupole to be consistent with zero at large scales when we 
remove the redshift space distortions. However, different trends are 
observed in the post-reconstruction mean correlation function, de¬ 
pending on the details of the reconstruction implementation and/or 
the simulations. The PTHALOS mocks show a slightly positive 
quadrupole post-reconstruction compared with QPM and PATCHY, 
which show a slightly negative quadrupole. We note that the imple¬ 
mentation used in PTHALOS is from [Padmanabhan et 2l.\{20\2) , 
while the one used in the QPM and PATCHY mocks is ours; the 
cosmology of the three different simulations is also different. A 
perfect reconstruction would remove large-scale anisotropy from 
the correlation function when we remove the redshift space distor¬ 
tions from the catalogues. However, in the case plotted, there is a 
residual anisotropy in the quadrupole at large scales (negative or 
positive), showing that the reconstruction is not perfect. The aim 
of this article is to disentangle the relation of this residual with the 
smoothing length and check the effect of the smoothing length in 
the BAO anisotropic post-reconstruction results. 

In this work, the main metric we use to evaluate the perfor¬ 
mance of reconstruction is the anisotropic BAO fits. Additionally, 
we explore the smoothing effects on the correlation functions and 
on the displacement field as a way of understanding the anisotropic 
fit results. We find that the smoothing length affects the quadrupole 
amplitude. Furthermore, we find that the differences in the ampli¬ 
tude also depend on the implementation. However, we show that 
the differences in the amplitude of the quadrupole do not deter¬ 
mine the anisotropic results. We also explore whether the origin of 
the improvement in anisotropic fits was related to a better estimate 
of the displacement field. 

This study addresses a crucial point in current BAO analysis, 
especially in the context of the final data release of BOSS galaxy 
data. The results we found are BOSS specific; however, we include 
a section with some additional tests to explore how these conclu¬ 
sions scale with the bias and number density that the results could 
be generalised to other surveys. 

The layout of this paper is as follows. We introduce the re¬ 
construction implementation in Section [^ We present the simu¬ 
lations used in our study in Section [^ and the analysis methods 
for anisotropic BAO measurements in Section [^ We then present 
the results in Sections [^ and Section [^presents the effect of the 
smoothing in anisotropic BAO fitting results and Sectionj^presents 
the effect of smoothing on the displacement estimation accuracy. 
We conclude in Section [7] 


2 BASIC RECONSTRUCTION ALGORITHM 

The algorithm of density field reconstruction has been described in 
Eisenstein et al.[([2007[), [Padmanabhan et al.[([2012[) and [Burden et[ 
ai( |20T4) . We describe the most general algorithm applied to hi- 
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Table 1. Smoothing scales tested with simulations 


Reference 





R’s (/i-i Mpc) 

Best R 

observable 

Mock catalogues 

Padmanabhan et al. 

(2012 

1 


10,15,20, 25 

15 

So(r) 

N-body 

Burden et al. 

mu 

; Burd 

en et al 

i 2015| 

5,8 10,15 20 ,40 

15 

P{k) 

2LPT Mocks 

Burden et al 

^015} 




5 10,15 

- 


2LPT Mocks 

Noh, White 1 

Padmanabhan (2009 

5,10, 15, 20, 25, 30 

- 

^(r) 

N-body 



Figure 1. Performance of reconstruction tested in different kinds of sky- 
type mock catalogues. We show the mean of monopole [top panel], 
quadrupole [bottom panel] from 100 mocks pre-reconstruction and post¬ 
reconstruction. The different colours are different kinds of mock cata¬ 
logues: darker shades are pre-reconstruction and lighter shades are post¬ 
reconstruction catalogues. In blue are t he PTHALOS (Manera et al.|2013) . 
In green are the Quick Particle Mesh jWhite, Tinker & McBride|2014t . 
In red are the PATCHY (Kitaura et al 2015; companion paper). For the 
monopole we get pretty similar results in the BAO fitting range. However, in 
all cases the post-reconstruction quadrupole is not exactly zero; there is an 
extra correlation (negative or positive). The purpose of this work is to disen¬ 
tangle the relation of this residual with the smoothing length and check the 
effect of the smoothing length on the BAO anisotropic post-reconstruction 
results. 


ased tracers considering redshift space distortions, the angular and 
radial mask of the data.We focus on the part where the smoothing 
scale enters the reconstruction algorithm. The algorithm proposed 
by |Eisenstein et al.| ( |2007| ) can be directly summarised as follows: 

• Estimate the over-density field from the galaxy positions us¬ 
ing an interpolation method. We are using the Nearest Grid Point 
(NGP) interpolation method. 

• Smooth the over-density field using a Gaussian filter with 
smoothing scale R in order to eliminate high k non-linearities 
(small scales). 

WG(fc) =exp(-i?2fcV2) (1) 

• Solve the Eq. ^ using the Zeldovich approximation. 

+ /V-(^-7^)7^ = --^, (2) 

where ^ is the displacement field, f the galaxy position, Sg the 
galaxy contrast of density and b the bias. 

We can solve Eq. ^ in configuration space following a finite 


differences approach ( [Padmanabhan et al.|2012| l or in Fourier space 
( [Burden et al.||2014] |. In order to solve it in Fourier space, it is 
assumed that ■ r)r 

is irrotational, which enables us to approximate it as the gradient 
of a potential field, i.e solve Eq. The implications of this ap¬ 
proximation were explored in [Burden et aI| ( [20T5] ) showing that the 
displacement field is underestimated when the irrotational compo¬ 
nent is neglected, suggesting an empirical formula to correct this 
effect. 

In our implementation we also use the Fourier transform method 
to solve for the displacements. However, we follow a simpler ap¬ 
proach by neglecting the effect of the RSD when measuring the 
density field, leading to this equation : 

v-'5 = -^, (3) 

b 

instead of solving equation]^ , and we estimate the displacement 
field by: 


^ = IFFT 


—ikSg{k) 

¥b 


(4) 


We verify that this choice does not affect the conclusions of the 
tests performed (see Appendix [A]). The effect of applying different 
RSD corrections on the displacement field is performed in [Burden] 
[et al[ ( [2015] ). A study of the effects of applying different RSD cor¬ 
rections in terms of anisotropic fits is performed in Vargas-Magana 
(in preparation). 

• Once the displacement is computed, we move the particles 
positions by the corresponding displacement, — $, to approximate 
their initial Lagrangian positions. This step provides the short-scale 
modes of the reconstructed density jPad manabhan, White & Cohn| 
[2009|[Noh, White & Padmanabhan[200"^ 

• Move an additional "^rsd we want to eliminate the redshift 
space distortions at large scales in the catalogue: 

^RSD = -/(^ • r) f. (5) 


• Generate a uniform random sample and move the particles us¬ 
ing the displacement field previously estimated from data. This step 
provides us the large-scale modes of the reconstructed density. 


The 2-point correlation post-reconstruction is then defined as: 

RR{r, ii) 

where D accounts for the data, S for the “shifted” random sample 
and R for the “non-shifted” random sample. The DD describes the 
pair counts per r-/i bin with data-data, SS the same for the shifted 
random pair-counts and RR for the non-shifted randoms. DS will 
be the pair-counts per r-/i bin taking one point from data and one 
from the shifted random set. 

The current reconstruction algorithm treats separately the 
small and large scale modes. While small scale modes stay in the 
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shifted galaxy field denoted by ”D”, the large scale modes are im¬ 
printed in the shifted random catalogue, denoted by ”S”. The re¬ 
constructed density field is then defined as 6rec = Sd — Ss that 
represents the sum of the two contributions of the large- and small- 
scale modes ( [Padmanabhan, White & Cohn|2009| l. 


3 SIMULATIONS 

In this work we use two kinds of mock catalogues. One set of ap¬ 
proximate mock catalogues that enables us to have a sufficiently 
large number of realizations to test the BAO anisotropic fitting re¬ 
sults. A second set, composed of a small number of high fidelity 
mocks based on N-body simulations, guarantees that the velocities 
are more accurate, enabling us to test the accuracy of the recon¬ 
structed velocity field for cosmological applications. 


3.1 Quick Particle Mesh Mocks. 

Quick Particle Mesh (QPM hereafter) mocks were generated for 
BOSS clustering analysis. These mock catalogues use low mass 
and force resolution particle-mesh simulations employing 1280^ 
particles in a (2560/i“^Mpc)^ box run with large time steps. At 
select times, the particles and their local density smoothed on 
2/i“^Mpc scales were dumped; these particles were then sampled 
(with a density-dependent probability) to form a set of mock halos 
that are then populated using a halo occupation distribution ( [White, | 
[Tinker & McBride|2014| ). 

We use “Sky mocks,” which match the observed number den¬ 
sity of BOSS galaxies and follow the radial and angular selection 
functions. These mocks are required to study the anisotropic galaxy 
clustering and to extract conclusions applicable to current analy¬ 
sis of BOSS-DR12 galaxy samples. We used the version of QPM 
mocks that matches CMASS North Galactic Cap samples for Data 
Release 12 of BOSS. 


by a and the anisotropic warping of the coordinates parametrized 
by e. a and e parametrize the geometrical distortion derived from 
assuming a “wrong” cosmology when calculating the galaxy corre¬ 
lation function. 

The parameters a and e are defined as 



2/3 1/3 

Cf. — ol ^ ^11 ^ 

(7) 

and 


(8) 


\aAj 

where a± and a\\ 

are defined by 



r_L = Q(_Lr_L,obs 

(9) 


^11 ~ Q(|©||,obs- 

(10) 


Here, the ‘obs” subscript denotes the observed coordinates. The 
coordinates without ”obs” subscripts correspond to the assumed 
cosmology. r± and r\\ are respectively the transverse and parallel 
to the line line-of-sight galaxy separations. 

The transverse shift, a±, allows us to measure DA{z)/rs, 
where Da{z) is the angular diameter distance to redshift and rs 
is the sound horizon scale. The line-of-sight shift, allows us to 
measure cz/{H{z)rs), where H{z) is the Hubble parameter. This 
is done using 


and 


_ DA(z)rs,fid 
D^Vs 


_ H^'^(z)rs,fid 
“ Hiz)rs 


( 11 ) 


( 12 ) 


Note that when analyzing the mocks in the cosmology in which 
they are define, a = 1 and e = 0. 


3.2 RunPB Simulations 

In Sectionj^ we use RunPB Simulations (RunPB hereafter) ( [White[ 
[et al.[[20T5] ) for the study of the velocity field. The catalogues 
are based on high-resolution realisations of the A—CDM model 
with Qrn = 0.292 and h = 0.69, employing 2048^ particles 
in a periodic box of side length 1380/i“^Mpc for a total volume 
of 2.6/i“^Gpc^. The values used at redshift z = 0.55 for the 
growth factor and the power spectrum amplitude are / = 0.76 
and (78 = 0.62. The simulations were run with the TreePM code; 
the mock catalogues are described further in [White[ ( |2010[ ). Briefiy, 
halos were found using the friends-of-friends algorithm. We use a 
cut in the halo mass {Mhaio > 1.10^^ M©) which mimics the 
bias property of the CMASS survey.We can do that under the as¬ 
sumption that most of the CMASS galaxies are Brightest Central 
Galaxies (BCG) (90%) and so have the same velocities as the ha¬ 
los. 


4 METHODOLOGY 
4.1 Anisotropic Analysis 

In this paper, we follow the multipole fitting procedure described 
in [Xu et aLl ( [2012[ ) and [Anderson et al.| ( [2013[ ), which extracts mea¬ 
surements of the isotropic dilation of the coordinates parametrized 


4.2 Nonlinear Models for the Correlation Function 

The template for the 2D nonlinear power spectrum, following from 
[Fisher et al| ( [1994[ ) reads as follows: 

P{k, /i) = (1 + /i, J:s)PNL(k, ii). (13) 

where the term (1 + /3/i^) corresponds to the Kaiser 
model for large-scale redshifts distortions, which produce an 
anisotropic damping and PnL{ k) is the nonlinear power spectrum. 
F{k, /i, Ss) is the streaming model for FoG given by: 

F{k,fl,T,s) = ;, 2 ^ 25 ] 2 ) ’ 

where is the streaming scale. In this work, we consider the de- 
wiggled power spectrum Pdw To get the templates for the mul¬ 
tipoles, we decompose the 2D power spectrum into its Legendre 
moments: 

2 / +1 

PiAk) = J ^ Pt{k,ij)Li{ix)dix, (15) 

which can then be transformed to configuration space using 

= J ^ Pt,di{kr), (16) 

where ji{kr) is the 1-th spherical Bessel function and Li(/d) is the 
1-th Legendre polynomial. 
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4.3 De-wiggled Template 

The de-wiggled template is a power spectrum prescription widely 
used in clustering analysis ( |Xu et aL|2012[ [Anderson et al.|20T^ 
|2013| ). This phenomenological prescription takes a linear power 
spectrum template to which we add the nonlinear growth of struc¬ 
ture. The de-wiggled power spectrum is defined as: 


Pdw{k,fd) = [Plin{k) — Pnw{k)] 


xexp 




+ Pn 


(17) 


where Piin{k) is the linear theory power spectrum and Pnw{k) is 
a power spectrum without the acoustic oscillations ([Eisenstein &| 


of the standard Gaussian damping of BAO, T>nl ■ 

Snl = (Sfi + Ei)/2. (18) 


Hu (|1998|)). E|| and are the radial and transverse components 


E atl models the degradation of signal because of nonlinear growth 
of structure. 


4.3.1 Multipole Fitting 

In order to measure a and e, we define a fiducial cosmology and 
then fit template multipoles (monopole and quadrupole) of the cor¬ 
relation function in this fiducial cosmology to the observed data. 
We use the model described in |Xu et ar] ( |2012| ) to generate the tem¬ 
plates for multipoles of the correlation function. This model is sum¬ 
marized in detai l in|Xu et al.|(|2012|),|Vargas-Magana et al.| ( [2014| ), 
[Anderson et al.| ( [^i2| ) and |Anderson et al.| ( |2013| ). The models we 

fit to our observed multipoles, ^o{r) and ^ 2 (^), are: 

fo(r) = Bo^o.obsC?-) + Ao{r), 

6(»’) = ^.obsW + A2{r), (19) 


where 


+ = ( 20 ) 

These A^{r) terms are used to marginalize out broadband (shape) 
information through the ... a ^,3 nuisance parameters. We use 
a fit he monopole and the quadrupole in a range of 50 < r < 
200/i“^Mpc, with 40 bins of 8 Mpc each. Knowing our model 
uses 10 we have finally 30 degrees of freedom for the fit. 

In order to find the best-fitting values for a and e, we minimize 
the function 


X^ = {rn-dfC ^{m-d), (21) 


where rh is the model vector and d is the vector of data. We scale 
the inverse sample covariance matrix, C~^, using 


^—1 _ ^—1 -Smocks Abins 2 

" A^ocks - 1 ■ 


( 22 ) 


to correct t he bias in the estim ate of the true inverse covariance 
matrix C~^ { [Hartlap et al.|2007] l. 

Error estimates for a and e are obtained by walking a grid in 
these two parameters to map out the likelihood surface. Assuming 
that the likelihood surface is Gaussian, this allows us to estimate 
a a and ae as the standard deviations of the marginalized ID likeli¬ 
hoods of a and e respectively. 


Table 2. Interpolation parameters 


Simulation 

Interpolation 

Method 

Box Size 
(fe-iMpc) 

Grid 

Pix Size 
(fe-iMpc) 

QPM Sky 

NGP 

3400 

512 

6.6 

RunPB Mocks 

NGP 

1380 

512 

2.7 


5 SMOOTHING EFFECT ON ANISOTROPIC BAO 
ANALYSIS. 

First, we study the effect of the smoothing scale in the BAO analy¬ 
sis described in section|^ The results presented in this section were 
obtained using “sky mocks” as described in section [3.1[ We ap¬ 
ply the reconstruction algorithm described in section [fusing four 
different smoothing scales: R=5, 10, 15 and 40/i~^Mpc. Then, we 
compute the correlation function of the reconstructed catalogues. 
Finally we apply the BAO fitting methodology and evaluate recon¬ 
struction performance using the fits on a and e parameters and their 
respective uncertainties as a metric. I^We summarize the details of 
the interpolation method used in the density estimation in Table 


5.1 Multipoles Results 

In Figure we show the mean multipoles from 200 QPM re¬ 
constructed mocks with the four smoothing scales. The upper 
panel is for the mean monopole and the lower panel for the mean 
quadrupole. The shaded regions correspond to the square root of 
the diagonal elements of the covariance matrix. 

For the monopole we observe two effects: 1) At large scales 
we find, as expected, that the different smoothing scales affect 
the level of sharpening of the mean monopole. The cases 5, 10, 
and 15 Mpc/h show an enhanced peak compared with the pre¬ 
reconstruction monopole. However the R=5 is slightly lower com¬ 
pared to the 10-15 smoothing scales. The R=40 is clearly decreas¬ 
ing the contrast in the BAO peak. 2) At small scales we find that the 
different smoothing scales affect the shape of the monopole. These 
scales are not important for the fitting of the BAO feature. The dif¬ 
ferences at those scales are mostly related to the residual redshift 
space distortions. 

Concerning the quadrupole, a perfect reconstruction would 
show a mean quadrupole consistent with zero (without any large- 
scale anisotropies). Figure illustrates that the typical smooth¬ 
ing scale of 15 h~^ Mpc gives a quadrupole not completely con¬ 
sistent with zero. The disagreement appears worse when using a 
large smoothing scale where the correlation increases in the inter¬ 
val 50 — 100/i“^Mpc. Using a smaller smoothing scale reduces the 
negative correlation observed in the quadrupole to be in agreement 
with zero within the fitting range. 

In summary, just by comparing the results in terms of correla¬ 
tion functions, it is not clear that the reconstruction is performing 


^ The effective smoothing scale is a convolution of the Gaussian smoothing 
and the grid smoothing. Hereafter we refer only to the Gaussian smoothing. 
We considered that the grid smoothing is not significantly changing the ef¬ 
fective smoothing. In the ID case, using a Gaussian smoothing equal to a 
grid size, the effective smoothing is ~ 15% larger for NGP, or ~ 30% 
larger for CIC. The 5 Mpc smoothing, assuming a 5 Mpc grid, is really ~ 6 
Mpc for NGP and 6.5 Mpc for CIC, which is not very different. 
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better with any smoothing scale. Two effects affect differently the 
BAO anisotropic fits: 

(i) The precision of a is driven by the sharpening of the BAO 
feature. Thus we expect, that 10-15, gives best precision on a fol¬ 
lowed very closely by 5, and that 40 smoothing scale increases the 
error on a. 

(ii) The precision on e depends on the capability that our nui¬ 
sance parameters (Bo, Ao, 2 ) to absorb the residual quadrupole. We 
can expect that a quadrupole consistent with zero is likely to give a 
good fit, thus we expect that 5 Mpc/h will provide best constraints 
on e, followed by 10 and 15. 

Our next task is to test these hypotheses and check how the smooth¬ 
ing scales affect the anisotropic fits. 

5.2 Anisotropic Fits on the Mean Correlation Function 

In order to determine which smoothing scale is performing better, 
we fit the nonlinear damping parameters and the streaming 
parameter using the mean information over the mocks. In Ta¬ 
ble [3] we summarize the results for different values of R. We find a 
minimum value for Ey^^ for a smoothing scale of 5 h~^Mpc, in¬ 
dicating that the mean correlation function is less nonlinear using 
this filter. When we increase the smoothing scale, the value Ey^^ 
increases for 10 and 15 Mpc/h, indicating that the function is be¬ 
coming more non linear. We note that the values we get for the 
Ey^^ for the 40 case are smaller compared to 15. At 40 Mpc/h we 
expect to have a poor performance of reconstruction, since we lose 
information as we are smoothing scales that are in linear regime. 
We suspect that the unexpected values obtained for the nonlinear 
damping are generated because we are using the fitting template 
for the post-reconstruction case, i.e. we are assuming that the non¬ 
linear damping is isotropic. But this is most likely not a good ap¬ 
proximation, since with this smoothing we are not removing most 
of the nonlinear evolution of the density. Thus, it would be more 
accurate to do the fitting with both damping parameters free in¬ 
stead . In order to test that the assumption Ey = E^ is break¬ 
ing down for R=40, we performed for all cases a fit without this 
assumption (see bottom panel of Table [^. Following this method¬ 
ology we found that the get the more symmetric values for R=5, 
however the smaller value for Eatl is given by 10 Mpc/h followed 
by R=5 Mpc/h. And the larger value is coming for R=40 Mpc/h. 
We notice also that the best fits gives very asymmetric values. Thus 
probably indicating that post-reconstruction is not necessarily well 
described by this assumption. The lower values of just indicates 
that the data is too good to be true, thus usually can indicate two 
different possibilities: 1) our model is valid but that a statistically 
improbable excursion of 2) we overestimate the errors. We do 
not consider our errors overestimated as they are coming from the 
sample variance of the simulations renormalized by the \/Nsim- 
We also notice that the /3 value fitted is the smallest for 
5/i“^Mpc, indicating that the correction for the Kaiser effect is less 
important for a smaller smoothing scale. In the case of the stream¬ 
ing parameter Eg, which is related to the Fingers-of-God (FOG) ef¬ 
fect, we do not expect to find any improvement (to find any change) 
in the best fitting value for this parameter, since reconstruction does 
not take care of virial velocities. 

5.3 Anisotropic Fits Results 

In Figure]^ we show the results from fitting 200 QPM mocks re¬ 
constructed using the four smoothing scales. We show the mean 


Table 3. Values for the fixed parameters [Sy = E^ 
(/i“^Mpc),Es(/i“^Mpc) , Bo and /3] fitted to the mean of 200 
QPM mocks varying the smoothing scale of reconstruction R (/i“^Mpc) 


R 

E|| = Ex 

E. 

X^ld.o.f 

Bo 

/9 


5 

4.9 

0.0 

18.4/30 

0.98 

-0.06 


10 

5.7 

0.0 

9.4/30 

1.01 

-0.12 


15 

6.5 

0.0 

14.3/30 

1.28 

-0.11 


40 

5.9 

-6.0 

11.3/30 

1.31 

0.29 


R 

sii 

E± 

Eatl 

X^/d.o.f 

Bo 

/3 

5 

4.8 

5.6 

5.2 

1.0/30 

0.9 

0.01 

10 

5.9 

1.5 

4.3 

8.2/30 

1.0 

-0.13 

15 

6.9 

3.6 

5.5 

10.3/30 

1.2 

-0.10 

40 

5.6 

00 

bo 

7.3 

14.7/30 

1.3 

0.35 


value from the best fits for BAO-related parameters a [top panel], 
e [bottom panel]. We also show in Figure [8] results for ^ [bottom 
panel], and the [top panel] results for the fits. The dotted line 
indicates the expected values from the mocks. The quantitative fit¬ 
ting results are summarized in Tablej^for the best fits in anisotropic 
BAO analysis. 

The accuracy in the a value is very similar between 5 and 
15/i“^Mpc, but visibly degrades when going to higher smoothing 
scales. The smaller dispersion for a is for 10-15 h~^Mpc; however, 
the 15 /i“^Mpc is significantly biased (5.3 cr), thus, 10/i“^Mpc 
is the best option if we considered only a. As we find a relative 
large bias in 15 h~^Mpch~^Mpc compared to 5,10,40/i“^Mpc, 
we perform a cross-check between the results using the 4 smooth¬ 
ing scales as a sanity check. We included in Figurej^the dispersions 
plots of the q;i 5 fits from the mocks, and the R=5,10, and 40 Mpc/h. 
The legend includes the values of the correlation coefficient in the 
three cases. The results show that effectively we are using the same 
set of mocks as they show very hight correlation coefficient and 
that changing the smoothing scale of reconstruction affects slightly 
the correlation between the results. 

For e, the least significant bias (bejae lower) is for the 5 
h~^Mpc smoothing scale. Considering a and e simultaneously, 
the best option is 5 h~^Mpc which gives the less significant bias 
with smaller dispersion (0.8cr in a and 0.9cr in e ). 

We find maximal variations in the mean value for best fit pa¬ 
rameters, ~0.5% for a and 0.3% for e, producing AD a, AH ~ 
0.5%. In the case of /3, the precision seems to be very similar for 
the 5 and 15/i“^Mpc smoothing scales and slightly lower for the 40 
h~^Mpc. The accuracy is also similar for 5 and 15h~^ smoothing 
scales and clearly degrades when using a higher smoothing scale. 
The Id.o.f values are very similar between the different smooth¬ 
ing scales explorec0 

Figure [^presents the distributions of the uncertainties for a 
and e parameters for different smoothing scales and the quantitative 
fitting results for the uncertainties are summarised in Table We 
can see that the smoothing scale strongly affects the error distribu¬ 
tions. and that the error bars decrease using 5 h~^Mpc smoothing 
scale. The mean error is reduced by 13% for a and 24% for e when 
passing from 15 to 5 h~^Mpc. 

^ The values of the x^ are similar to those obtained with the recon¬ 
struction implementation used in previous analysis. We add more discussion 
about the results using different implementations in AppendixlAl 
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Figure 2. Mean of 200 QPM mocks NGC reconstructed with different smoothing scale 5h“^Mpc [green], 10h“^Mpc [cyan], 15h“^Mpc [magenta] and 
40/i“^Mpc [blue]. The smoothing scale is correlated with the negative correlation observed in the quadrupole; a smaller smoothing scale decreases the 
correlation observed; a 5h“^Mpc smoothing scale erases the quadrupole almost completely. Right panel: Zoom to the monopole [top panel] and quadrupole 
[bottom panel] in the BAO range. 


1.10 


1.05 


g 1.00 


0.95 
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Figure 3. Dispersion plots of 0:15 vs q; 5 ,io, 40 - The legend includes the 
values of the correlation coefficient in the three cases. The large values of 
the correlation coefficient just shows that we are effectively using the same 
set of mocks for the test and that changing the smoothing scale of recon¬ 
struction affects the correlation between the results slightly, but remains 
large. 


. 5 MPc/h 0.756922 
. lOMpc/h 0.939044 
. 40Mpc/h 0.790126 



• • 


In this section, we have shown how the smoothing scale affects 
the anisotropic clustering results and we found that the smaller bias 
and error bars are obtained using a smoothing scale of 5/i“^Mpc. In 
order to associate a systematic error, we used the scale 15/i“^Mpc 
as the fiducial smoothing. We observe that varying the smoothing 
length from 15 to 5/i“^Mpc generates a Aa = 0.005 and Ae = 
0.001 (Table 1^. Concerning the errors, the variation observed is 



Figure 4. Mean of best fits for ol and e for 200 QPM mocks NGC analysed 
with different smoothing scales. The error bars are given by the standard 
deviation from 200 realisations. For a, the smaller dispersion is for 10-15 
h~^Mpc; however, the 15 Mpc/h is significantly biased (53 a). The best 
is 5 h“^Mpc, which has less significant bias with small dispersion (O.Scr). 
For e, the less significant bias (belcTe lower) is for the 5 h“^Mpc smoothing 
scale given their small bias and dispersion. 


Aaa = 0.002 and Acre = 0.012. Expressing these variations in 
the best fits and their uncertainties to the final measurements of the 
angular diameter distance and Hubble parameter, we get variations 
of Acfda ^ 40% and Acfh ~33 %. 


5.4 Dependence on Reconstruction Implementation 

In Appendix we show that the results obtained in this work 
for the best smoothing scale, in terms of anisotropic fits perfor¬ 
mance, are independent of the particular implementation of recon- 
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Figure 5. Mean of best fits for /3 and ^ jd.o.f. for 200 QPM mocks 
analysed with different smoothing scales. The error bars are given by the 
standard deviation from 200 realisations. The 5,10 h~^ Mpc smoothing 
scale gives similar RMS and bias for linear redshift distortion parameter (3 
but slightly more bias for larger smoothing. Even though /3 is a nuisance 
parameter in BAO analysis, it is interesting to observe the value fitted, as 
it indicates the level at which the redshift corrections are using the right 
value of the velocity field. The /d.o.f values are very similar between 
the different smoothing scales explored. 


Table 4. Best fits for a and e. Mean and RMS from 200 reconstructed QPM 
mocks. Second block refers to the results using the covariance matrix from 
1000 mocks with a fixed smoothing scale of 15 Mpc/h. 


R 

a 

RMS 

bpi 

^ Cy. 

e 

RMS 

O-e 

/3 

RMS 

5 

0.9992 

0.0136 

0.8 

0.0007 

0.0109 

0.9 

-0.035 

0.090 

10 

1.0004 

0.0128 

0.5 

0.0024 

0.0142 

2.4 

-0.025 

0.088 

15 

1.0048 

0.0128 

5.3 

0.0021 

0.0169 

1.7 

-0.019 

0.085 

40 

1.0002 

0.0147 

0.2 

-0.0048 

0.0224 

3.0 

0.1020 

0.077 

5 

1.0003 

0.0138 

0.2 

0.0002 

0.0115 

0.5 

-0.028 

0.079 

10 

1.0012 

0.0137 

1.2 

0.0021 

0.0143 

2.1 

-0.042 

0.089 

15 

1.0027 

0.0142 

2.7 

0.0029 

0.0185 

2.2 

-0.030 

0.093 

40 

0.9980 

0.0158 

0.5 

-0.0011 

0.0256 

0.6 

0.1493 

0.102 


Table 5. Uncertainties in a and e parameters. Median, 16 and 84 percentiles 
of the uncertainties distributions from 200 reconstructed QPM mocks. 


R 

^oc 

A(7q:% 


b 

<1 

5 

0 01 Q«+0 0021 

U.UiC)D_Q QQJ^g 

-9.9 

n m 1 o-l-0.0028 

0 . 0 iiS_Q 

-45.9 

10 

0 0X42+®-°®^® 

U.Ui^Z_Q QQ^g 

-5.9 

U.UiO/ _Q QQ 24 

-23.4 

15 

0 01 ^1 +0 0030 

U.UiOi_Q QQ 23 

- 

n 0910+0.0062 

U.Uzi«_Q QQ3g 

- 

40 

0 01 c; 0 +O OO 42 

U.UiDU_Q QQ 23 

-0.7 

n n 9 QQ+ 0.0080 
U.UZOO_Q QQ 43 ^ 

6.9 
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Figure 6. Histograms of uncertainties in a and e measured for the indi¬ 
vidual reconstructed mocks for different smoothing scales, 5 (green), 10 
(cyan), 15 (red), and 40 (blue). The distributions depend strongly on the 
smoothing scale used in the reconstruction. 


struction algorithm. We compare our implementation of the recon¬ 
struction algorithm to the method implemented by |Padmanabhan| 
et al.|{|201^. We compare results in terms of the multipoles. Figure 


Al shows the mean multipoles for the different implementations 
using three different smoothing lengths, 5, 10, and 15 h~^ Mpc. 
We observe that the monopole behaviour is consistent between dif¬ 
ferent implementations in the fitting range. The differences ob¬ 
served are only important at scales smaller than 20 h~^ Mpc. 
The amplitude of the quadrupole between different implementa¬ 
tions is significantly different. RP implementation shows a positive 
quadrupole which means that the reconstruction implementation is 
over-correcting the anisotropy. On the other hand, our implementa¬ 
tion of reconstruction shows a negative quadrupole, which means 
that our implementation is under-correcting the anisotropy; how¬ 
ever, both quadrupoles are very similar in shape. Differences in the 
quadrupole are generated by two effects, redshift distortions cor¬ 
rection and effects of the angular and radial selection functions, 
which are implemented in slightly different ways. Further explo¬ 
ration and quantification of these two contributions which gener¬ 
ate differences in the quadrupole amplitude are treated in Vargas- 
Magana et al (in prep). 

We fitted the multipoles post-reconstruction from both imple¬ 
mentations and compared the results. The results are completely 
consistent between both implementations, the differences between 
the mean of the two implementations is 0.1% for a and e. The RMS 
of the best fits are similar at 0.1% in both quantities. Concerning 
the errors, the dispersion of the distributions are also very similar; 
however, there is a systematic larger error in RP implementation 
compared to RV in a and e. However, the same trends are observed 
in both implementations; the error in e is monotonically decreasing 
as we apply a smaller smoothing scale. In the case of a the errors 
we get with the 5 and 10 h~^ Mpc smoothing scale are very simi¬ 
lar, but smaller than using the 15 Mpc that is regularly applied 
in BAO analysis. 


5.5 Dependence on Covariance Noise 

We also test the effect of the covariance noise on the fitting results. 
The fitting results presented in previous sections use the covariance 
matrix generated from the 200 reconstructed mocks with a differ¬ 
ent smoothing scale. In this section, we substitute this covariance 
matrix by the covariance generated with 1000 mocks with the 15 
Mpc/h smoothing scale and we fit the fourth sets of multipoles post¬ 
reconstruction of 200 with different smoothing scales. This choice 
is motivated by previous results [Vargas-Magana et al.| ( |2016| ) indi- 
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Figure 7. Mean of best fits for a and e for 200 QPM mocks NGC analysed 
with different smoothing scales fitted with the 200-covariance matrix using 
the reconstructed mocks with different smoothing scales [in red] compared 
with the fits performed with the 1000-covariance matrix generated from 
post-reconstrucion mocks using a fixed smoothing scale of 15h“^Mpc [in 
blue]. The error bars are given by the standard deviation from the 200 re¬ 
alizations. For a, the smaller dispersion is for 5-10 h“^Mpc. The best is 5 
h“^Mpc, which has less significant bias with small dispersion (0.2cr). For 
e, the less significant bias (be/fJe lower) is for the 5 /i“^Mpc smoothing 
scale given their small bias and dispersion. 



Figure 8. Mean of best fits for [3 and jd.o.f. for 200 QPM mocks 
analysed with different smoothing scales. The error bars are given by the 
standard deviation from 200 realisations. The 5,10 h~^ Mpc smoothing 
scale gives similar RMS and bias for linear redshift distortion parameter (3 
but slightly more bias for larger smoothing. Even though /3 is a nuisance 
parameter in BAO analysis, it is interesting to observe the value fitted, as 
it indicates the level at which the redshift corrections are using the right 
value of the velocity field. The "x^/d.o.f values are very similar between 
the different smoothing scales explored. 


eating the covariance noise could affect the stability of the fitting 
results. By performing this change we are testing the impact in the 
fitting results just from the noise observed in the covariance. Even 
this procedure is neglecting the effect of smoothing scale in the co- 
variance post-reconstruction, we consider it to be testing the confi¬ 
dence we can have in the fitting results given the reduced number 
of mocks we used for these tests. 

Results are shown in the second block of Table |4] and also 
in Figure [7] The results follow similar trends than using the 
“noisy” covariance matrix (even though the fourth significant fig¬ 
ure changes slightly). The less biased results for a and e are ob¬ 
tained from the smaller smoothing scale 5 Mpe/h followed by the 
10 Mpe/h and 15 Mpe/h smoothing scales and then degrades for 40 
Mpe/h. The dispersions shows variations for some cases by ~0.001 
related the error bars themselves are noisy estimates, i.e. with 200 
mocks there should be more scatter in the error bars than there 
would be with 1,000 mocks. 

An interesting feature of the fits performed with the 1000- 
covariance is that the significant bias observed in the 15 Mpe/h case 
for the previous sections reduces from 5.3cr to just 1.1 g which in¬ 
dicates the noise in the covariance was generating this large bias. 
We highlight that even though the variations derived from the noise 
in the covariance, the main conclusions about the performance of 
the fits for the different smoothing scales remain unchanged. The 
smaller smoothing scale of 5 Mpe/h is giving the best results. Sum¬ 
marizing the results, the maximal variations in the mean value 
for best fit parameters considering this new covariance reduces to 
0.002 for a and e, producing AD a, AH ~ 0.003 and 0.004%. 


6 SMOOTHING EFFECT ON ACCURACY OF 
DISPLACEMENT FIELD ESTIMATION. 

In this section, we explore the effect of the smoothing scale on the 
reconstructed displacement field to figure out the origin of the im¬ 
provement observed in the post-reconstruction anisotropic cluster¬ 
ing. In particular, we investigate if the improvement is related to 
the better estimation of the displacement field when we reduce the 
smoothing scale applied in the reconstruction algorithm. 

Moreover, the peculiar velocity of galaxies is a valuable quan¬ 
tity in cosmology, since it contains complementary information to 
that enclosed in the galaxies’ positions. In the literature, there are 
many strategies with different approximations to obtain the velocity 
field. Reconstruction provides the simplest approach to get a model 
dependent reliable velocity field at large scales. For reconstruction, 
the velocity field is obtained using the well-known Zeldovich ap¬ 
proximation, in which the displacement field is given by the gra¬ 
dient of the gravitational potential at q. The range of applicability 
of this method is limited to very large scales, as it fails to describe 
the dynamics of a nonlinear field ( [Kitaura & Angulo|2012| l. Even 
though it has a limited range of applicability, the velocity field from 
reconstruction provides a direct method to constrain gravitational 
model comparing with velocity measurements 

In this section, we analyze the accuracy of the velocity field 
derived from reconstruction and then the effect of the smoothing 
length on the estimated velocity field. Instead of working with a 
velocity field we transformed the velocity to displacements using 
the continuity equation: 

= (23) 

where / is the growth factor at redshift of the snapshot and H{z) 
is the Hubble parameter. 

First we evaluate the quality of the recovered displacement 
field from reconstruction compared with simulations. This test is 
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important for validating the use of displacement field for other cos¬ 
mological applications such as the measurement of the kSZ effect 


(|Schaan et al.|2015[ Planck Collaboration, Ade, P. A, R,, Aghanim,| 

|N., et al.|2015||Hernandez-Monteagudo et al.|2015| l. Secondly, we 

present an estimate of the uncertainties associated with the recon¬ 
structed displacements for individual galaxies and the bias of the 
displacement estimations with respect to the true displacements as 
a function of the smoothing length. 

For this section we use a different set of mock catalogues: the 
cubic RunPB real space mocks. We use RunPB mocks because they 
are expected to have more accurate velocities; even if we do not 
have a large number of these mocks to study the statistical prop¬ 
erties for BAO fits, we have enough galaxies in each of them to 
accurately compare the individual velocities with the reconstructed 
fielcQ 


6.1 Reconstructed Displacement Field 

The purpose of this subsection is to show that reconstruction is 
providing a reliable estimate of the true displacement field but it 
is limited to reproduce the large scales. We analyze the quality of 
the reconstructed displacement field as compared to the true dis¬ 
placement field. We compare the individual displacement of each 
galaxies with the reconstructed one in the correspondi^ position 
in a grid of 276^ pixels of size (5 h~^ Mpc)^. In Figure]^ we show 
the x-component of the displacement field, from reconstruction 

(right panels) and the true corresponding displacement from simu¬ 
lations (left panels) using a Gaussian filter of 10 /i“^Mpc. We only 
show one component of the vector field, as the other two compo¬ 
nents are quite similar. The is a 3D scalar field; for the illustra¬ 
tion, we show the average value in a slice of 50 h~^Mpc over one 
direction. The plot shows the reconstructed displacement field re¬ 
produces the simulation displacement field at scales of 10/i“^Mpc, 
demonstrating the accuracy of the reconstructed displacement field. 
However, we observe that, at small scales, the reconstructed dis¬ 
placement field shows fewer structures compared with the mock 
catalogue displacement field. 


6.2 Error and Bias of Reconstructed Displacements 

Up to now, we have been using the values of the displacement field 
in pixels. Now we are interested in comparing the individual values 
of the displacement for each galaxy to evaluate their agreement to 
the simulations. We expect only be correct on the linear component 
of the velocity field. We are also interested in associating an error 
with those velocities; this is interesting if we use those velocities 
for inferring cosmological information for applications other than 
BAO. 

We study the dispersion between true and reconstructed dis¬ 
placements, for which we estimate the 2D histograms mock-by¬ 
mock for different smoothing scales using the RunPB cubic mocks 
(Figure [T^. 

The dotted black line indicates the 1-1 relation. The solid line 
comes from a 2D Gaussian fitting to the 2D histogram; it represents 
the angle between the major axis and the ordinate axis. This angle 
cannot be interpreted as the bias of the reconstructed displacements 


Table 6. RunPB Cubic Error and bias from equations and[^ from 
1 reconstructed RunPB Cubic mock with different smoothing lengths. The 
best results are obtained for the 5h~^ Mpc smoothing scale. 


R 


a 

Gn 

OL 

5 

0.87 

1.01 

2.23 

2.21 

10 

0.83 

0.77 

2.08 

2.70 

15 

0.77 

0.64 

2.05 

3.20 

40 

0.60 

0.33 

1.74 

5.27 


compared with the true displacements, but it provides an illustra¬ 
tion of the effect. In Section [F.2.1| we describe the methodology we 
follow to characterise the bias and noise in the displacement esti¬ 
mation and we present the results in |6.2.2| The derivation of the 
equations is presented in Appendix [B] 


6.2.1 Theoretical Estimation of Bias and Noise 

In order to get a sensible estimate of the bias and noise level, we 
follow an approach in terms of probability distributions. We write 
the reconstructed displacement as a biased estimation of the ’’true” 
displacement ^ plus a Gaussian noise: 

= a X + A/'(0, ctat), (24) 

where a is the bias term and aN characterises entirely the noise. 
Because the noise is uncorrelated to the displacement, we then 
have: 


~2 2 2,2 

aR — a cr^+CTiv, (25) 

where o-r is the variance of the “reconstructed” displacements for 
the smoothing scale R and cr^ is the variance of the “true” displace¬ 
ments. We must notice we do not have direct access to a and aN- 
The observable quantities are the variance or standard deviation 
from and 'fR (i.e. a^ and gr) and a the angle of the correlation 
between and 

While the equation ( [^ provides a relation between the ob¬ 
servable quantities and the intrinsic bias and noise, we still need to 
break the degeneracy between the two contributions. We propose 
using the scatter plot between '^r and f). Following the calculus of 
Appendix we get the following expressions for the bias of the 
reconstructed displacements and the associated noise (a, gn) for a 
smoothing scale R given the measurable quantities (a, g^, 


~ 9 ? 

aa^ 


(26) 


2 

Gn 


~2 


^CTr- 




(27) 


The dispersion of the de-biased displacement is then given by 
Gn/ oi. This quantity enables us to evaluate the quality of the esti¬ 
mator depending only on the smoothing scale R, gn, and a. 


^ The ideal set of mocks for testing the displacement field accuracy would 
be high fidelity mocks with reliable velocities with the properties of the 
survey; however, we only have available high fidelity mocks without the 
mask of the simulations (RunPB mocks). 


^ gr using directly the standard deviation of the reconstructed values. In 
the case of simulations, we have access to which is obviously not the 
case using data. 














Effect of smoothing of density field on Reconstruction. 


11 


Simulation 



Ho.25 

H 0.20 

- 0.15 

- 0.10 

- 0.05 

- 0.00 

- -0.05 

-0.10 



Reconstruction 



Figure 9. Comparison of the displacement field from a RunPB simulation [left] and from reconstruction [right]. The reconstructed field reproduces the 
structure of the true velocity field, demonstrating the accuracy of the reconstructed displacement field. Both of the fields are smoothed with a 10 Mpc/h 
Gaussian filter. The color bars are in Mpc/h. 


R=5Mpc/h 


R=10Mpc/h 


R=15Mpc/h 


R=40Mpc/h 


CO 






Figure 10. RunPB mocks. Scatter plots between true and reconstructed displacements from different smoothing scales, from left to right R=5,10,15,40 h~^ 
Mpc. The mocks used in this test were in real space. The dashed white line indicates the 1-1 relation. The solid line comes from fitting a 2D-Gaussian to the 
2D histogram; it represents the angle between the major axis and the ordinate axis. This angle does not indicate the bias of the reconstructed displacements 
compared with the true displacements but provides an illustration of the effect. The best result is obtained for the 5 h~^ Mpc according to our quantitative 
results shown in Table[^ These mocks have more realistic velocities compared to QPM mocks. 


6.2.2 Measurement of Bias and Noise 

We measure the cr^^ and cr^ standard deviation from reconstructed 
and simulation displacements. We measure the correlation coeffi¬ 
cient 'ijj). We calculate the biased coefficient a as: 


(TffCifijdR.'f) 

We deduce the two quantities, the bias a and the noise, gn, using 
equations and Once we have the values for the bias and 
the noise, we can determine which smoothing scale is the best. We 
expect a result similar to the correlation coefficient with extra in¬ 
formation from the contribution of the intrinsic noise, cr at, for the 
different smoothing scales. We expect a lower noise for the larger 
smoothing scale, as we are deleting the nonlinear scales more ef¬ 
ficiently. However, we also expect to have a poorer precision in 
the displacement reconstruction because we are losing some struc¬ 
tures corresponding to the linear theory. Therefore, we expect a 
lower value for the bias a as we increase the smoothing scale. The 
best smoothing scale is therefore the one that balances the two ef¬ 


fects and corresponds to the minimum value of the ratio cta^/q; as 
demonstrated in the appendix. 

The results are shown in Table for RunPB cubic mocks in 
real space. Figure summarises bias and noise results. The top 
panel shows the correlation coefficient < > as a func¬ 

tion of the smoothing scale, the intermediate panel shows the qual¬ 
ity factor gn/ oL as a function of the smoothing scale. Finally, the 
bottom panel shows the noise g and the bias a as a function of the 
smoothing scale. 

The best results are obtained for the highest correlation coef¬ 
ficient and the lowest quality factor. We find that the quality value 
is a monotonic function of the smoothing scale. The ratio crAr/a is 
minimized for R — Mpc and also that the maximal correla¬ 
tion coefficient value {C = 0.87) corresponds to the R = 

Mpc. We can see also from the bottom panel of Figure[^ that both 
the bias, a, and the noise, gn, are decreasing with the smoothing 
scale, but at a different ratio. As shown in the appendix B, the ratio 
is the way to measure the accuracy of the estimator (equation B.IO). 
FigurepTjalso allows to explain why R — Mpc is performing 
better, while a smaller smoothing results in a improved correlation. 
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Figure 11. Top panel shows the correlation coefficient < > 

, central panel the quality factor ajy/a, and bottom panel the bias a and 
the noise a, as a function of the smoothing scale. The best results for cubic 
mocks RunPB are for 5 h~^ Mpc, while the best result for the Sky mocks 
are obtained for 10-15 h~^ Mpc. 

the noise level increases more abruptly than the correlation, thus 
giving a poorer displacement field. From the top and intermediate 
panels of Figure [TT] we can see that the correlation coefficient and 
the ratio cta^/q; are symmetric and so provide equivalent informa¬ 
tion. However, while the correlation coefficient provides a single 
information, with the methodology we proposed, we have access 
to the bias and the noise of the estimator separately. This point is 
important when one wants to use the reconstruction methodology 
to derive the velocity field directly and associate errors to it. 

6.2.3 Halo Bias Dependence 

In this section, we explore how the conclusions we got for BOSS- 
like samples in the last section scale for other biases and number 
densities so the results could be applied to other surveys. We show 
the effects in terms of the displacement field as the bias effect in 
the correlation function; anisotropic fits are treated in another ap¬ 
pendix. For this test, we generate several halo samples from RunPB 
simulation, applying different halo mass cuts. The sample informa¬ 
tion is summarised in Table we show the number density, shot 
noise level, and halo bias. The halo bias is determined by compar¬ 
ing the power spectrum of halos with the dark matter power spec¬ 
trum from the simulations. 

We show the scatter plots between true and reconstructed 
displacements for RunPB Mocks in Figure Each line of 
plots represents a different halo cut, i.e represents a sample 
with a different halo bias. From top to bottom, Mhaio = 
[lel2, 5el2,9el2, lel3, 2el3,4el3, 5el3]/i"^M0. The different 
columns show the different smoothing scales, from left to right 
R=[2,5,10,15,20,40]/i“^ Mpc. The dashed black line indicates the 
1-1 relation. The solid line comes from fitting a 2D Gaussian dis¬ 
tribution to the 2D histogram; it represents the angle between the 
major axis and the ordinate axis. This angle does not indicate the 


Table 7. RunPB mocks samples generated from different halo mass cuts. 


Mcut 

number density 

Shot Noise 

Halo Bias 

kmax 

1.0el2 

0.00287320 

348.044 

1.29069 

0.568604 

3.0el2 

0.00105167 

950.872 

1.55206 

0.364857 

5.0el2 

0.000616345 

1622.47 

1.71838 

0.326531 

7.0el2 

0.000425321 

2351.17 

1.85922 

0.261644 

9.0el2 

0.000319666 

3128.27 

1.96803 

0.234149 

l.OelS 

0.000281704 

3549.83 

2.02379 

0.234149 

2.0el3 

0.000117975 

8476.39 

2.41207 

0.167890 

3.0el3 

6.69662e-05 

14932.9 

2.71752 

0.134476 

4.0el3 

4.35019e-05 

22987.5 

2.95412 

0.107754 

5.0el3 

3.05197e-05 

32765.7 

3.14878 

0.0964309 


bias of the reconstructed displacements compared with the true dis¬ 
placements but provides an illustration of the effect. 

The results of the bias and noise (equations[^and|^ are pre¬ 
sented in Table The quality factor and coefficient of correlation 
indicate what is the best smoothing scale for each case. The best re¬ 
sults are for the higher correlation coefficient and the lower quality 
factor. Figure summarises bias and noise results for the sam¬ 
ples generated with different halo mass cuts. The left panel shows 
the correlation coefficient < C{^r^ 'ip) > and the right panel the 
quality factor crA^/a as a function of the smoothing scale. 

The results for the correlation coefficient and the quality fac¬ 
tor show two regimes: 1) for a small smoothing scale the result 
depends strongly on the bias; 2) for a large smoothing scale the re¬ 
sults converge in the correlation coefficient and quality factor for 
all the biases. These two behaviours could be understood consid¬ 
ering number density and shot noise. The greater the mass cut is, 
the greater the bias is and also the greater the shot noise level of 
the samples. The samples that are more biased are more affected 
by the shot noise; for those cases, a smoothing scale from 5-10/i“^ 
Mpc is preferable, while for low-biased tracers, a smaller smooth¬ 
ing scale, from 2-5h~^ Mpc, also gives a good result. On the other 
hand we can see that the result is similar for all the biases for the 
large smoothing scale, indicating that the smoothing is compensat¬ 
ing for the shot noise effect on the reconstructed displacements. 


7 CONCLUSION 

We summarize the results presented in this paper. 

• We test the reconstruction algorithm using four different 
smoothing scales, 5,10,15, 40 Mpc/h with QPM BOSS-like mock 
catalogues and we study the effect in BAO anisotropic analysis. 
We find that the different smoothing scales affect the multipoles. 
We observe variations in the amplitude of the monopole at low 
scales and in the sharpening of the BAO feature. The effects on 
the quadrupole are larger than on the monopole. The smoothing 
scale of 15 h~^ Mpc produces a quadrupole not completely con¬ 
sistent with zero. Taking a smaller smoothing scale reduces the 
negative correlation observed in the quadrupole to be in agree¬ 
ment with zero within the fitting range. A large smoothing scale 
increases the negative correlation in the quadrupole at scales be¬ 
tween 50 — 100/i“^Mpc. 

• We show that the smoothing scale affects the anisotropic clus¬ 
tering results. The results indicate that the best choice for the 
smoothing length, in terms of anisotropic clustering, is given by the 
smaller smoothing scale of Mpc. This smoothing scale shows 
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R=2Mpc/h R=5Mpc/h R=10Mpc/h R=15Mpc/h R=20Mpc/h R=40Mpc/h 





w Reconstruction (Mpc/h) 
R=2Mpc/h 


w Reconstruction (Mpc/h) 
R=5Mpc/h 


w Reconstruction (Mpc/h) 
R=10Mpc/h 


w Reconstruction (Mpc/h) 
R=15Mpc/h 


w Reconstruction (Mpc/h) 
R=20Mpc/h 


w Reconstruction (Mpc/h) 
R=40Mpc/h 




Figure 12. We show the scatter plots between true and reconstructed displacements for RunPB mocks. Each line of plots represents different halo cuts (i.e 
represents a sample with a different halo bias). The different columns represent the different smoothing scales, from left to right R=2, 5, 10, 15, 20, 40h~^ 
Mpc. The mocks used in this test were in real space. The dashed black line indicates the 1-1 relation. The solid line comes from fitting a 2D Gaussian 
distribution to the 2D histogram, representing the angle between the major axis and the ordinate axis. This angle does not indicate the bias of the reconstructed 
displacements compared with the true displacements but provides an illustration of the effect. The results of the bias and noise are presented in the Table 
The higher correlation coefficient and the lower quality factor indicates the best result for each case. 
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Figure 13. The left panel shows the correlation coefficient < C{'ipR, 'ip) > and the right panel the quality factor (Jiv/Q; as a function of the smoothing scale. 
The different samples follow similar trends: the best results are obtained for the higher correlation coefficient and the lower quality factor; using the 5 h~^ 
Mpc smoothing scale gives the best results for the most of the cases. The more biased samples are more affected by the shot noise, and for those cases, a 
smoothing scale from 5-10h~^ Mpc is preferable, while for low-biased tracers, a smaller smoothing scale, from 2-5h~^ Mpc, also gives a good result. On 
the other hand, we can see that the result is similar for all the biases for the large smoothing scale, indicating that the smoothing is compensating for the shot 
noise effect on the reconstructed displacements. 


the smaller bias and error on a an e error bars simultaneously. The 
variations in the mean value for best fits parameters using 15h~^ 
Mpc compared with Mpc are ~0.5% for a and 0.3 for e, 
when considering the noise in the covariance the variations in a 
are 0.3 and 0.2 for e. The variations produce AD a ~ 0.3 and for 
AH ~ 0.4%. The smoothing scale affects the precision at which 
we can measure the BAO parameters. Taking 15/i“^Mpc as a ref¬ 
erence, the median error using 5 Mpc/h is reduced by 10% for a 
and 45% for e. The effect on the uncertainties of Da{z) and H{z) 
are 40% and 30% respectively. The variations in the uncertainties 
does not consider the covariance noise. 

• We explore the effect of the smoothing scale in the recon¬ 
structed displacement field to investigate our hypothesis, whether 
the improvement observed in the post-reconstruction anisotropic 
clustering was coming from a better estimation of the displacement 
field when we reduce the smoothing scale. We estimated the cor¬ 
relation coefficient between the two fields and also we provided 
a methodology to estimate the noise in addition to the correlation 
coefficient, through the quality factor (defined as the ratio gn/ol). 
The correlation coefficient and the quality factor are symmetric and 
so provide equivalent information. However, while the correlation 
coefficient provides a single information, with the methodology we 
proposed, we have access to the bias and the noise of the estimator 
separately. This point is important when one wants to use the recon¬ 
struction methodology to derive the velocity field directly and as¬ 
sociate individual errors to it. We find that the quality value, on /ol 
is minimized for R — Mpc and also that the maximal corre¬ 
lation coefficient value {C — 0.87) corresponds to the R — 

Mpc. 

• We explore how the conclusions we get for BOSS-like sam¬ 
ples scale for other bias and number densities so the results could 
be applied to other surveys. We show the effects in terms of the 
displacement field (and in the appendix, the bias effect in the cor¬ 
relation function). For this test, we generate several halo samples 
from RunPB simulation, applying different halo mass cuts. The re¬ 
sults for the correlation coefficient and the quality factor show two 
regimes: 1) For small smoothing scales, the result depends strongly 


on the bias; the samples that are more biased are more affected by 
the shot noise; for those cases, a smoothing scale from 5-10/i“^ 
Mpc is preferable, while for low-biased tracers, a smaller smooth¬ 
ing scale, from 2-5h~^ Mpc also gives a good result. 2) For large 
smoothing scales, the results converge in the correlation coefficient 
and quality factor for all the biases, indicating that the smoothing 
is compensating for the shot noise effect on the reconstructed dis¬ 
placements, and the dependence on bias is weaker. 

Further work should be done to provide a theoretical framework for 
the empirical results shown in this article. The optimal smoothing 
scale could vary for different tracers (different kinds of galaxies) 
or could possibly also depend on the redshift and environmental 
properties. A future work will focus on exploring these different 
dependencies to look for ways to improve the reconstruction tech¬ 
nique for future surveys such as cBOSS, DESI, EUCLID. 
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Table 8. RunPB mocks samples for several halo mass cuts. Error and bias 
from reconstructed RunPB cubic mocks with different smoothing lengths. 


R 


a 

O'N 

cx 

Mass cut = 1.0el2 

5 

0.87 

0.87 

1.94 

2.21 

10 

0.81 

0.70 

1.95 

2.79 

15 

0.76 

0.58 

1.93 

3.31 

20 

0.72 

0.49 

1.88 

3.79 

40 

0.59 

0.29 

1.59 

5.41 

Mass cut = 3.0el2 

5 

0.87 

0.91 

2.00 

2.19 

10 

0.82 

0.72 

1.97 

2.74 

15 

0.77 

0.59 

1.95 

Zll 

20 

0.72 

0.50 

1.90 

3.75 

40 

0.59 

0.29 

1.61 

5.38 

Mass cut = 5.0el2 

5 

0.87 

0.94 

2.06 

2.18 

10 

0.82 

0.73 

2.00 

2Jl 

15 

0.77 

0.60 

1.96 

3.23 

20 

0.73 

0.51 

1.91 

3.70 

40 

0.59 

0.30 

1.61 

5.31 

Mass cut = 7.0el2 

5 

0.87 

0.97 

2.13 

2.19 

10 

0.82 

0.75 

2.03 

2.71 

15 

0.77 

0.61 

2.00 

3.23 

20 

0.72 

0.52 

1.94 

3.71 

40 

0.59 

0.30 

1.65 

5.35 

Mass cut = 9.0el2 

5 

0.87 

0.99 

2.19 

2.20 

10 

0.82 

0.76 

2.06 

2.69 

15 

0.77 

0.62 

2.02 

3.22 

20 

0.72 

0.53 

1.96 

3.70 

40 

0.59 

0.31 

1.67 

5.33 
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APPENDIX A: COMPARISON OF DIFFERENT 
RECONSTRUCTION IMPLEMENTATIONS. 

In this section we show that the results obtained in this work are 
independent of the particular implementation of the reconstruction 
algorithm. We compare our implementation of the reconstruction 
algorithm (hereafter “RV”) to the method implemented by |Pad-| 
[manabhan et al.| ( |2012| ), hereafter denoted “RP.” The differences 
between the two methods could be summarised as follows: 

• Solution of Displacement Field. RV implementation solves the 
Poisson equation in Fourier Space assuming a real space density 
field, RP in configuration space using a finite differences technique. 

• Redshift Space Corrections. RV neglects the redshift space 
corrections, while RP takes them into account. This choice is tested 
to give the accuracy required for BAO anisotropic analysis at sub¬ 
percent precision in Vargas-Magana (in prep). The corrections are 
not important because the effect on the quadrupole is a change of 
the amplitude, and the peak position and contrast are not affected. 

• Survey Mask and Density Estimate. Both methods follow dif¬ 
ferent approaches for dealing with the effect of geometry. 

We used the same 100 mocks described before and we run recon¬ 
struction with RP and RV implementation using the same param¬ 
eters. We use for this test a bias=2.1 instead of bias=1.87 because 
the reconstructed catalogues were available only for this value of 
the bias for the RP implementations. We compare results in terms 
of the multipoles and anisotropic fits. We show in Figure [aT] the 
mean multipoles for the different implementations using three dif¬ 
ferent smoothing lengths, 5, 10, and 15 h~^ Mpc. We observe that 
the monopole behaviour is pretty similar between different imple¬ 
mentations in the range of the fitting; the differences observed are 
only important at scales smaller than 20 h~^ Mpc. In the case 
of the quadrupole, the amplitude of the quadrupole between dif¬ 
ferent implementations is significantly different. The quadrupole 
for RP implementation seems to over-correct the anisotropy, while 
it seems to be under-correcting in our implementation; however, 
both quadrupoles are very similar in shape. Differences in the 




r (Mpc/h) r (Mpc/h) 

Figure Al. Multipoles for different smoothing scales: 5,10 and 15 h~^ 
Mpc for different implementations. In our implementation: 5 h~^ Mpc in 
cyan, 10 h~^ Mpc in blue and 15 h~^ Mpc in black. For Padmanabhan et 
al.’s implementation: 5 h~^ Mpc in yellow, 10 h~^ Mpc in green and 15 
h~^ Mpc in red. In the monopole, we observe that, in the fitting region, both 
implementations perform similarly; however, there are some differences at 
very small scales (< 20h~^ Mpc). For the quadrupole, we observe impor¬ 
tant differences in the amplitude. The quadrupole for Padmanabhan et al. 
seems to be kind of over-correcting, while in our implementation seems to 
be under-correcting; however, both quadrupoles are very similar in shape. 
Differences in the quadrupole are generated by two effects: redshift dis¬ 
tortions correction and the effects of the angular and radial selection func¬ 
tions that are implemented in slightly different ways. Further exploration 
and quantification of these two contributions, which generate differences in 
the quadrupole amplitude, are treated in Vargas-Magana et al. (in prep). 

quadrupole are generated by two effects: redshift distortions correc¬ 
tion and the effects of the angular and radial selection functions that 
are implemented in slightly different ways. Further exploration and 
quantification of these two contributions, which generate differ¬ 
ences in the quadrupole amplitude are addresses in Vargas-Magana 
et al. (in prep). 

Tables [AT] and |A2| summarise the best fitting results. The first 
table refers to the best fits of a and e as well as their variance. In the 
second, table we show the mean values of the errors and the RMS 
of the error distributions. Additionally, Figures [A^ and | A3 [ show the 
comparison between both implementations’ best fit parameters and 
errors (mean and RMS). The differences between the mean of two 
different implementations is ~ 0.1% for a and e, indicating the 
results are completely consistent between both implementations. 
The RMS of the best fits are also similar at 0.1% in both quantities. 

Concerning the errors, the dispersion of the distributions are 
about the same; however, it seems to be a systematic larger error 
in RP implementation, as compared to RV in a and e. Although 
the trends are observed in both implementations, the error in e is 
monotonically decreasing as we apply a smaller smoothing scale. 
In the case of a the errors we get with the 5 and 10 h~^ Mpc 
smoothing scales are very similar to but smaller than using the 15 
h~^ Mpc that is regularly applied in BAO analysis. 

Finally, in Figure [A4| we put together the results from Section 
|5.3| and the results from this section for the anisotropic fits using 
the same methodology but varying the bias from b=1.87 to b=2.1. 
The results show that the bias is not affecting the best fits, i.e, the 
conclusions we get for the smoothing scales in terms of anisotropic 
results are valid even considering other bias values. It is interesting 
to mention that even when comparing Figures[^and | A1 1 we observe 
that the bias affects the quadrupole amplitude post-reconstruction, 
this does not affect the best fitting values. The uncertainties are 
shown in Figure [A5] for both bias values. We observe that the dis¬ 
tributions of the uncertainties are quite similar for 10-15 h~^ Mpc 
smoothing scale for both a and e distributions; however, the lower 
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Table Al. Best fits. Mean and RMS from 100 reconstructed QPM NGC 
mocks for different Reconstruction Implementations. 


E 

a 

RMS 

e 

RMS 

This work 

5 

0.9961 

0.0091 

0.0036 

0.0098 

10 

0.9980 

0.0096 

0.0023 

0.0107 

15 

1.0040 

0.0095 

0.0050 

0.0135 

Padmanabhan et al 

5 

0.9970 

0.0100 

-0.0002 

0.0088 

10 

0.9991 

0.0116 

0.0022 

0.0112 

15 

1.0053 

0.0107 

0.0036 

0.0141 
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Figure A3. Smoothing Fitting Results.Right, mean era and RMS of a a 
distribution [right]. Left mean cfe and RMS for both implementations. We 
use a value of the bias=2.1 different than the value used in the article, i.e 
bias=1.87 


Table A2. Uncertainties on the best fits. Mean, RMS from 100 recon¬ 
structed QPM NGC mocks different Reconstruction Implementations. 


S 

<7’a 

RMS 

CTe 

RMS 

This work. 

5 

0.0158 

0.0054 

0.0194 

0.0082 

10 

0.0151 

0.0044 

0.0199 

0.0071 

15 

0.0163 

0.0051 

0.0271 

0.0084 

Padmanabhan et al 

5 

0.0171 

0.0063 

0.0186 

0.0090 

10 

0.0185 

0.0049 

0.0212 

0.0070 

15 

0.0186 

0.0051 

0.0277 

0.0076 
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Figure A4. Smoothing Fitting Results. Top panel, mean a and RMS of a 
distribution [right]. Bottom panel, mean e and RMS for two values of the 
bias=1.87 and bias=2.1. Best fits and dispersion are similar for both values 
of the bias for all the smoothing scales. 


smoothing scale shows a large dispersion using the bias=1.87 com¬ 
pared to the bias=2.1. This result seems to show that the bias can 
play a role in the uncertainties’ determination when using small 
smoothing scales. 


degeneracy between the bias a and the noise effect. We propose to 
use the correlation plot between and We will express the 


APPENDIX B: BIAS AND NOISE DETERMINATION 

While the equation ( [24| provides us a relation between the vari¬ 
ances of reconstructed and true displacement, we need to break the 
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Figure A2. Smoothing Fitting Results. Right panel, mean a and RMS of a 
distribution [right]. Left panel, mean e and RMS for both implementations. 
We use a value of the bias=2.1 different than the value used in the article, 
i.e bias=1.87 
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Figure A5. Smoothing Fitting Results. Top panel, mean cfa and RMS of 
a a distribution [right]. Bottom panel, mean erg and RMS cfe distribution 
for for two values of the bias=1.87 and bias=2.1. There is a large difference 
in the dispersion of the a and e distributions for the lower smoothing scale; 
however, for 10-15 h~^ Mpe, the results are consistent between both bias 
values. 
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theoretical probability distribution knowin^the bias a, 

the noise standard deviation ctat. 


v{i)R,i)) = v{tp) X (Bi) 

Assuming a Gaussian distribution, we can express the probability 
distribution of the “true” displacements by: 

Using the relation expressed in equation ( [24| we can get the 
expression for the conditional probability distribution V{'ipR, ^): 




1 


‘2i 7V(Tip (T TV 


exp 


{^R - a'lpy 


^ N 


(B3) 

The relation we expect to measure between 'ijjR and ^ is where 
V{'ipR^ 'ip) is maximum depending on ^r value, i.e. we are inter¬ 
ested in the value i^max, which maximises the exponential term for 
a given ^jr. Then we use the partial derivation: 


d I (tpR - aipY 

+ a% 


= 0 . 


(B4) 


'^R 


Developing this calculus, we obtain the following expression 


for 'Iprr 


and so a is given by: 


'4^max — 


a'lpR 

t2 ’ 

^1P 


(B5) 


^R = a'll^max ^ a = - 

a 


(B6) 


We find the relation between ^jr and 'ijjmaxi which can be 
thought of as representing the apparently most correlated value. 

The 2D correlation provides us an estimation of the bias a 
which is different than a except for a null noise term ctat = 0. 

We have two independent equations (Eq. and Eq. |B6| ) 
which connect (a, gn) to the measurable quantities (a, cr^, cr^^). 

Combining Eq.|B6|and Eq.l^ we obtain: 


2 ’ 


(B7) 


and 


2 _ ~2 CFr 

— CTR — ~9 2 


(B8) 


for whi ch, we present the comparison with measurements in Fig- 
The plot shows the Vi^ipR, 'ip) for a = 0.5, ctat = 8 and 


Bl 


Gip = 10. The dashed black line represents the relation between 
^R and ^ without noise term (so a), the solid black line repre¬ 
sents the theoretical relation obtained using Eq. |B6| (so a); the solid 
red line represent the major axis orientation using a 2D-Gaussian 
fit. The yellow diamonds correspond respectively to the measured 
'ipmaxi^R) on thc left and the mean value (^{^r) ) on the right 



Figure Bl. V{'ipR, 'ip) for a = 0.5, = 8 and cr-^ = 10. The dashed 

black line represents the relation between 'ipR and pj without noise term (so 
a), the solid black line represents the theoretical relation obtained using Eq. 
|B6| (so d); the solid red line represents the major axis orientation using a 
2D-Gaussian fit. The yellow diamonds correspond respectively to the mea¬ 
sured 'ipmaxi'ipR) on the left and the mean value (^'ip{'ipR)'^ on the right 
panel. 



Figure B2. V{'ipR,'ip) for alpha = 0.8, gn = 8 and = 10. The 
dashed black line represents the relation between ^r and 'ip without noise 
term (so a), the solid black line represents the theoretical relation obtained 
using Eq. |B6| (so d); the solid red line represents the major axis orientation 
using a 2D-Gaussian fit. The yellow diamonds correspond respectively to 
the measured 'ipmax{^ r) on the left and the mean value (yji'ipR)^ on the 
right panel. 


panel. From these tests, we conclude that the bias and noise estima¬ 
tion described below describes the bias and noise levels observed 
in our Monte-Carlo realisations. 


We can write the de-biased estimator 'ipR inverting Eq. (241: 
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A/'(0,crAT) 

a 


(B9) 


Also we can write the dispersion on the de-biased reconstruction 
as: 


a 


(BIO) 


Gip is a measured quantity and so is known too. 

























